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Abstract. The quasistationary spreading of a circular liquid drop on a solid 
substrate typically obeys the so-called Tanner law, with the instantaneous base 
radius R{t) growing with time as i? ~ t^^^'^ - an effect of the dominant role 
of capillary forces for a small-sized droplet. However, for droplets of nematic 
liquid crystals, a faster spreading law sets in at long times, so that R t" with 
a significantly larger than the Tanner exponent 1/10. In the framework of the 
thin film model (or lubrication approximation), we describe this "acceleration" 
as a transition to a qualitatively different spreading regime driven by a strong 
substrate-liquid interaction specific to nematics (antagonistic anchoring at the 
interfaces). The numerical solution of the thin film equation agrees well with 
the available experimental data for nematics, even though the non-Newtonian 
rheology has yet to be taken into account. Thus we complement the theory of 
spreading with a post-Tanner stage, noting that the spreading process can be 
expected to cross over from the usual capillarity-dominated stage to a regime 
where the whole reservoir becomes a diffusive film in the sense of Derjaguin. 
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1. Introduction 

The spreading of liquid drops and films on a solid surface can be described by universal, 
"macroscopic" laws [Ij as soon as the thickness of the drop or film exceeds a few 
tens of nanometers. One such law is the so-called Tanner law, characteristic of the 
spontaneous spreading of small non- volatile drops on a flat substrate in a situation of 
complete wetting (see flgure[T|). After an initial transient regime, the base radius R 
of such a drop grows as i? ~ The law has been derived analytically [llISlll] and 

confirmed experimentally on many accounts [HI [H [U . The fundamental argument is 
that the hydrodynamics in the bulk of a drop are driven by capillary forces alone, which 
directly yields R ~ assuming a self-similar shape for the bulk, in the lubrication 
approximation O |3j. Alternatively, the trend can be regarded as a competition 
between the hydrodynamic dissipation (primarily in the contact Hue region of the 
drop) and an unbalanced capillary force [H [H [HI . 

The Tanner law is quite robust and typically offers an accurate description of 
the life of a droplet - which spans a few hours for liquids with moderate surface 
tensions and viscosities - apart from initial and final transients. The initial transient 
corresponds, e.g., to the deposit of the droplet on the substrate, and lasts less than 
a second for regular Hquids. As for the final state of spreading, for non-volatile 
droplets it is either a molecular film or a fiat, bounded structure - a "mesoscopic" 
pancake [H [TOl [HI [12] - which may be more favorable energetically than a molecular 
film. Pancakes occur when short-range substrate interactions promote dewetting, 
even though the overall situation is that of complete wetting [13] : although not very 
common, such structures have been observed experimentally [H]. The existence of 
a limiting configuration, with a finite value for the base radius R, implies that the 
late-time spreading dynamics typically slow down with respect to the Tanner law. 

By contrast to the trend of arrested spreading, it has recently been observed that 
the Tanner stage R ~ i^/^" can be followed by a faster spreading law R ^ t", with 
a > 1/10. Specifically, for spontaneously spreading nematic Hquid crystals [TSl fT6] . 
the value of a was found to be nearly twice as large as the exponent aTannor = 0.1 
characterizing the Tanner law, with a — 0.2 |15] and a = 0.19 [16,]. A more thorough 
analysis of the data suggests that the acceleration does not stop at a = 0.2: values 
as high as a = 0.3 can be observed at the end of the experiment. This "accelerating" 
trend is apparently in confiict with the notion of a Tanner stage terminated by the 
onset of a molecular film or equilibrium pancake, and its physical origin has yet to be 
clarified. 

We have already attempted a qualitative explanation of this post-Tanner trend 
in a macroscopic framework [IT]. In the present paper, our goal is to account for the 
acceleration quantitatively and for this purpose we resort to the well-accepted thin film 
model (TFM). Analyzing the thin film equation (TFE) we see that at late spreading 
times the disjoining pressure dominates capillary effects. Then the local thickness 
h{r,t) (see figure 1) obeys a diffusion equation, i.e., the whole droplet effectively 
becomes a diffusive film in the sense of Derjaguin [13]; the corresponding base radius R 
grows as R ^ t^/"^ . The experimentally observed transition from Tanner's law to power 
laws R ^ t" with a « 0.2 [151 [IS] in fact seems to be part of a crossover to a much 
faster spreading law than expected previously. In order to validate this observation 
we integrate the TFE numerically, extract relevant observables and compare their 
evolution to the experimental data. We find the general trend illustrated by the 
"numerical spreading" to be in good agreement with the spreading observed for nematic 
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Figure 1. Cross-section of a circular droplet spreading in a situation of complete 
wetting (cartoon). See section [TTI for details. 



droplets in [151 IIS- However, our model has yet to take into account the typically 
non-Newtonian rheology of nematic liquid crystals. Our work in progress will address 
this effect in separate publications. The macroscopic interpretation of the acceleration 
in terms of a negative line tension is also provided elsewhere [17] . 

Our paper is organized as follows. In section 2 we give an overview of our system of 
interest and summarize previously attempted explanations for its abnormal spreading 
behavior. Section 3 features a brief derivation of the thin film equation and presents 
key analytical results that are directly relevant to the problem of nematic droplets. 
Section 4 focuses on key properties of our numerical spreading process (essentially a 
brute-force solution of the TFE) , which is then compared quantitatively with physical 
experimental data in section 5. We conclude in section 6, relating our results to 
alternative concepts and providing outlook into our future work. 

2. Overview of the problem 

2.1. Anatomy of a spreading droplet: macroscopic versus mesoscopic 

Figure 1 represents the quasistationary state of a spreading droplet, which is composed 
of: a "macroscopic" liquid drop where shear is small and viscous forces are balanced 
primarily by variations of capillary Laplace pressure; a "mesoscopic" part subject to 
large shear, where viscous forces are balanced primarily by variations of disjoining 
pressure; a "microscopic" region featuring molecular precursor layers and a dry 
substrate [iHlII^. The relative sizes of these regions are not up to scale: the main 
purpose of figure 1 is to clearly distinguish the apparent contact line (typically inferred 
from the infiection point of the profile h{r, t) at a given time t) from the "real" contact 
line, governed by microscopic phenomena. 

Most analyses assume that the macroscopic and mesoscopic scales are well- 
separated, i.e., that the bulk of the drop is much wider than the mesoscopic "foot". 
Thus the bulk is well-approximated by a thin spherical cap with base radius R{t), 
contact angle 0{t) <^ 1 and nearly constant volume 



i.e., it is practically in equilibrium at constant volume Vcap and instantaneous base 
radius R{t). In this approximation it is also customary to assimilate R and with 



cap — 




(1) 
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Figure 2. Overview of ttie "accelerating" trend for a single spreading experiment 
(cyanobiphenyl 5CB droplet on silicon wafer). The quantities observed are iJsonm 
(microscopic contact line; small dots), -Rsoonm (second interference fringe; big 
dots) and Savcragod (slope averaged over 25 interference fringes; crosses). At 
t > 200 s the plot reveals the existence and growth of a large mesoscopic "foot" at 
the edge of the droplet, comparable in size to the base radius (see figures 1 and 
4a). It is also clear that both observed radii significantly deviate from Tanner's 
law R ~ t^/i'' and undergo a transient "acceleration", past both the R ^ t^/^ and 
R ^ i^/^'' power laws. The volume V of the droplet was not measured directly 
but it can be estimated as V = 1.236 x 10~^ mm'^ (see sections 14.41 and [5)l . 

their respective estimates inferred from the inflection point of h{r,t), i?inflection and 
^inflection (see figure 1), which can be obtained through optical observation of the 
apparent contact line [TH [16]. More immediately, these same optical experiments 
yield the location of the microscopic contact line (the characteristic thickness of which 
is 30 nm) and, e.g., i?3oonm(0 such that /i(i?3oonm, = 300 nm (see figures 2 and 
4a) . These measurements indicate that the length of the mesoscopic film ahead of the 
apparent contact line grows to millimetric sizes and becomes comparable with R at 
the end of the experiment. At this point the assumption of well-separated scales is 
clearly broken, and the macroscopic volume T4ap is significantly lower than the total 
volume 

/•oo 

V = 2tt I h{r,t) r dr. (2) 
Jo 

2.2. Post-Tanner spreading laws: experimental evidence and tentative explanations 

By contrast with previous evidence of Tanner's law, recent experimental studies of 
spontaneous spreading of nematic liquid crystals on hydrophilic [15] or hydrophobic 
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[T6] substrates revealed, after a transient Tanner stage, a surprising "acceleration" 
(actually a spreading process that "slows down more slowly" than the Tanner law). 
Figure 2 plots the relevant observables for one spreading experiment. The base radius 
R and contact angle were reportedly inferred from the inflection point of the thickness 
proflle (see flgure 1), although in fact the main observables were i?30nm and 6'avoraged 
(see figure 4a) rather than i?inflection and 6'inflection- Initially [Hj it was noted that 
Tanner's law crossed over to i? ~ with a « 0.2. It was also realized (see figure 4b 
in [H]) that the Tanner relation O'^ ^ Ca (where Ca = is the capillary number) 
does not hold for late spreading times: for small 9 and Ca one has 9 ^ Ca°'^^. The 
latter relation, together with the volume conservation condition R^9 ~ V, is consistent 
with R ~ Similar results were reported for spreading on hydrophobic substrates 
(see figure 6 in [16]), with 9 ~ Ca"'^ and R - 

The conclusions from these results are as follows. On one hand [15], direct 
estimates of a through R{t) are consistent with estimates via 9 (Ca), which apparently 
validates the hypothesis of an approximating spherical cap of constant volume 
(although the experimental data [H] offers no direct evidence for this). On the other 
hand, a is found to be significantly larger than the Tanner exponent axanncr = 0.1. 
This signifies that some unknown factor, other than the surface tension, comes into 
play. Moreover, similar values of a were obtained for different kinds of substrates 
[121 H^, which suggests that the acceleration is a robust effect rather than an artefact, 
and an intrinsic feature of nematic droplets. One should also note that the experiment 
in [1^ does not fully capture the post-Tanner transient and that exponents as high as 
a « 0.3 are observed at late times (see figure 2). 

Several qualitative arguments come to mind, which may or may not explain 
the reported acceleration. First of all, nematic crystals are known to have a non- 
Newtonian, shear-thinning rheology [2Ql[2l]. Shear thinning affects the fiow pattern 
and the dissipation rates and thus modifies the spreading dynamics, but it is not clear 
a priori whether the actual dynamics will be faster or slower than Tanner's power 
laws. 

A detailed analysis of the contact line dynamics in the framework of the thin 
film model [22] shows that characteristic shear rates in the capillary wedge and in 
the mesoscopic precursor decrease as the contact line slows down, and thus for a non- 
Newtonian fiuid the effective viscosity will increase with time. This corresponds to a 
modified spreading law R t" with a < 1/10. Numerical experiments were carried 
out [23] and confirmed a < 1/10 for shear-thinning fiuids and a > 1/10 for shear- 
thickening fiuids. Thus the dominant effect from shear thinning is that the spreading 
is slower than predicted by Tanner's law, and we must seek another mechanism to 
explain the "acceleration" observed for 5CB droplets. 

Among other factors that could be responsible for an acceleration of spreading, 
we should also cite: a) slippage at the substrate [24]; b) somewhat counterintuitively, 
densely distributed roughness, which for small or zero contact angles causes wicking 
and "enhances" the property of complete wetting [25]. Unfortunately, we must discard 
both these effects as possible causes of the observed acceleration in the case of 
our nematic droplets. First, slip has been shown to cause a logarithmically small 
contribution to the macroscopic spreading laws [24]; the acceleration observed in 
[121 [H] looks qualitatively different from a minor effect due to slip. As for roughness, 
the nematic droplets may encounter some anchoring defects on hydrophiHc substrates 
[15] , but hydrophobic substrates in p6] are definitely free of either chemical defects or 
topographic roughness; thus the consistent acceleration observed in both cases is not 



Post- Tanner spreading of nematic droplets 



6 



likely to be related to the "superwetting" properties of rough substrates. 

Finally, a very tempting approach is to describe the (non volatile) system in 
terms of its total free energy. The 1985 review by de Gennes pLj has explained 
Tanner's law in terms of an effective driving force (derived from the instantaneous 
free energy). The work of the driving force is balanced by dissipation, primarily 
hydrodynamic dissipation in the macroscopic "wedge" and mesoscopic "foot" in the 
vicinity of the apparent contact line. The assumptions of this approach can be 
challenged by hypothesizing incomplete dissipation in the foot/precursor [15] or by 
introducing the concept of dynamic line tension, which contributes to the unbalanced 
Young force and plays a dominant role at long spreading times |17j . 

We note, however, that line tension as an equilibrium concept is quite subtle 
[261 [27] and its generalization to a quasistationary situation should not be taken lightly. 
It is also hard to derive a consistent set of correction terms for the hydrodynamic 
dissipation. More generally, the notion of a driving force acting on the edge of a 
macroscopic, capillary drop - this notion breaks down when the size of the mesoscopic 
region ("foot") becomes comparable to that of the bulk of the droplet, which is 
apparently the case during the reported acceleration (this is indicated, e.g., by the 
evolution of i?30nm and i?3oonm on figure 2). This prompts us to describe the spreading 
droplet in a framework that resolves mesoscopic regions and does not use macroscopic 
approximations - the thin film model. 

3. Thin film equation: presentation and analytical results 

The thin film model (TFM) - related to both the "interface displacement model" 
and the "lubrication approximation" - is a continuum representation of the spreading 
dynamics, suitable for the study of thin films. While it may not accurately describe 
the spreading dynamics at molecular film thicknesses (see section \4A\\ . it is believed 
to work quite well for mesoscopic thicknesses, i.e., above several tens of nanometers. 
As opposed to macroscopic frameworks, the TFM accurately resolves quantities that 
would otherwise remain empirical, e.g., functionals of the thickness profiles in the 
vicinity of the apparent contact line and in the mesoscopic "foot" of a droplet. Notably, 
Tanner's law was consistently derived in the framework of the TFM by Voinov [2], 
Tanner [3] and de Gennes [H H]. Numerous authors have since used the TFM 
to validate, refine, and generalize the features of advancing contact lines and, by 
extension, the spreading dynamics of droplets [28l [29l [30j . 

3.1. Thin film model as applied to nematic droplets 

At the core of the framework is the thin film equation (TFE) . The simplest expression 
of the TFE is for a Newtonian fiuid with no slip at the substrate, in the approximation 
of small thickness gradients. In the following we briefiy derive a TFE for nematic 
droplets. 

We consider a quasistationary film of heterogeneous thickness h (x, y) covering a 
homogeneous, fiat substrate. Assuming that local equilibrium is achieved for all {x, y) 
and that the lateral fiows in the film have negligible inertia (low Reynolds number), 
we can write the following energy functional: 




(3) 
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Here a and (tsl are the nominal surface energies of a free interface and of a solid- 
liquid interface, respectively, whereas $ (h) is an effective interface potential acting as 
a correction to a + asL due to the fact that h is finite. As for the excess energy due 
to the curvature of the free interface, it is well approximated with §// (Vft.)^ dxdy 
in the small-slope approximation |V/i| ^ 1. It is common to neglect hydrostatic 
contributions to ([3]) in situations of complete wetting [H |30] . 

Considering £ [h] under volume-preserving variations of h, the quasistationary 
internal pressure p {x, y) is found to be of the intuitive form 

p = -crA/i - n (/i) , (4) 

which is a combination of the typical "capillary" Laplace pressure and of the 
"disjoining" pressure 11 (/i) = —dh^ih). The lateral pressure gradient Vp is relaxed 
through a so-called Poiseuille flow 

j = -|vp, (5) 

assuming a constant viscosity 77 and no sHp at the substrate (for a derivation, see [Ij 
or [30]). Finally, the conservation equation dth = —V • j, together with |[4]) and (O, 
yields 



a/i = -V • <^ - 



377 



-(jAh-Il{h) 

" v ' 

pressure p 



> . 



(6) 



lateral current j 

Due to the different nature of the two contributions to the pressure p, it is appropriate 
to rewrite ^ as 



dth 



■n 



V • [D{h)\/h] 



(7) 



where 



D = -- — — 



(8) 



h^ dn 

3r/ dh 

is the effective diffusion coefficient introduced by Derjaguin [13]. The expression of D 
contributes to the second-order term of ^ and plays a major role in the spreading 
dynamics at mesoscopic thicknesses. In a continuum representation, it is expedient to 
approximate the disjoining pressure n(ft.), and hence D{h), with a dominant long-range 
contribution, while introducing a phenomenological boundary condition (effectively a 
cutoff) in the nanometric range. Previous studies [H [4l [30] addressed the case of 
n (/i) = -^Ah~^, a single power law accounting for the cumulated effect of non- 
retarded van der Waals interactions for a film of finite thickness h {A being the 
Hamaker constant). However, in the case of antagonistically anchored nematic liquid 
crystals [HIH], the dominant term is 



I{=-K5^h-\ 
2 



(9) 



K being the bend-splay elastic constant and 5 the angle by which the director rotates 
over the thickness of the film, i.e., the difference between the anchoring angles at 
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both interfaces. Note that we only take into account the elastic energy in the bulk of 
the nematic; the anchoring energies (surface terms) are taken to be constant, i.e., we 
assume sufficiently strong anchoring at both interfaces with respect to the thickness 
of the film. Thus in our case IHl) reduces to 



D = 



KS^ 
3r/ 



(10) 



Interestingly, the elastic interaction typical of antagonistically anchored nematics 
yields a purely diffusive film in the sense of Derjaguin. To the best of our knowledge, 
this remarkable feature of a thickness-independent diffusion coefficient has not been 
emphasized previously. 

When describing the spreading droplet as a whole (as opposed to assimilating 
the apparent contact fine to a quasistationary hydrodynamic wedge [1, 4i30j), it is 
appropriate to rewrite |(7]) in a rotationally invariant geometry, i.e., with h depending 
only on the distance r to the vertical axis of the droplet (see figure 1), 



dth 



- — -dr I \h^rdr 
rj r [^3 



-dr 
r 



rdrh) 



D-dr 
r 



[rdrh) 



(11) 



Here we have used the fact that I? is a constant for nematics. 



3.2. Asymptotic spreading behavior 

Before we use l|lip for quantitative predictions, we can make an important qualitative 
remark about the two Hmiting cases of the TFE. Looking at the two terms on the 
right-hand side of ^ or ifTTj) . for sufficiently tall droplets, the disjoining pressure 
IT = ^K5'^h~'^ is neghgible with respect to the Laplace pressure over a large part 
of the droplet (the droplet is well-approximated by a spherical cap as in section 12.11 
and the Laplace pressure is 2a9/R ~ 27r [h{r ~ 0)] /14ap)- In this approximation, the 
TFE is essentially fourth-order and has the form 



^dth = --dr\h'^rdr 



Jt 

a r 



-dr [rdrh) 
r 



(12) 



Looking for self-similar solutions of the form h (r, t) — t {t "r), with the scaling 
chosen so that V = 2t: J h{r,t) r dr remains constant, lfT2|) yields a = 1/10, i.e.. 
Tanner's law [H [3l [Ml [29] . 

It is also clear that at late stages of spreading, as the droplet becomes fiatter, 
the Laplace pressure will eventually be dominated by the disjoining pressure, and 
the TFE [cf. ([7]) and ifTT])] will be essentially second-order, governed by the specific 
liquid-substrate interactions: 

dth^-dr[rD{h)drh]. (13) 
r 



In the nematic case, D has the constant expression (|T0)) and an obvious self-similar 
solution of (flSl) is a Gaussian bell defined by 

with an arbitrary origin for time. The base radius and contact angle as inferred from 
the infiection point are R — VSDt and 9 — so that the volume conservation 

relationship is y = as opposed to ^ = jR'^9 for the idealized Tanner regime. 
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Thus we see that the spreading process must cross over from an initial spreading 
phase, consistent with the generic Tanner's law, to another regime, specific to 
antagonistically anchored nematic liquid crystals: we shall see in section 14.41 that 



the characteristic time of the crossover scales as T = i^^V^ / D^j ■ The late-time 
evolution of a droplet is expected to be diffusive, i.e., measurements of R{t) will yield 
an "acceleration" from R ~ to i? ~ t^^^. The remaining problem is to establish 
the characteristics of the crossover and to compare it to the physical experiment. 

4. Numerical integration of the TFE: preliminaries 

We integrated ifTTj) in the form of a numerical spreading process, taking snapshots 
of the solution h{r,t) at preset time intervals (figure [3K). The numerical values used 
were ct = 30 x lO^^^.j^^-i^ 77 = 30 x lO^^p^.g and K6'^ = 12 x lO^^^ J-m^^ (i.e., 
D = 2.67 X 10~^°m^ • s~^). As for the boundary conditions, below /i = lA we 
extrapolate h{r,t) as an exponential tail (similar to a "maximal film" [U [U \30\ ) 
and ensure that the total volume V = 2Tr J h{r,t) r dr is preserved: in figure [3l 
V = 2.36 X 10~^ mm^. As for the initial conditions, we start with a perfect parabolic 
cap, to which we add a moderately smooth foot to avoid a computationally heavy 
singularity. Integration is explicit in time, with a fixed-size grid for r and an adaptive 
time step. 

Figure [SJd plots the spreading versus time in terms of the base radius R and 
contact angle 9, as inferred from the infiection point of h{r, t) (see figure 1 and figure 
[3fc). The volume of the numerical droplet is similar to the physical experiment p^, 
as well as the characteristics c = 30 x 10"'^ N • m~^, 77 = 30 x 10^"^ Pa • s and 
KS"^ = 12 10~^^ J • m~^. The acceleration agrees qualitatively with figure 2 (cf. figure 
5 in [H]): the crossover from the Tanner phase to significantly faster regimes occurs at 
times of the order of a minute, for a droplet of the same "caliber" as in |T6] (i? ~ 1 mm 
at i = 1 s) ; a characteristic time of the crossover corresponds to the intersection of the 
asymptotes in figure 3b, at i ~ 700s (about 12 minutes). 

Before comparing quantitatively the physical spreading experiment with our 
numerical resolution of the TFE, we shall voice a few words of caution about the 
applicability of the TFM. We shall also review the physical observables available to 
us. Finally, we shall discuss the issue of time origin in spreading experiments and the 
scaHng feature of equation (fTTj) . 

4.1. TFM applicability 

The oblique crosses on figures[3K and[3l: indicate the points at which the solution h (r, t) 
falls below e = 30 nm, which is the characteristic thickness at the edge of the physical 
mesoscopic precursor [TB]. For t = 2000 s (dash-dotted profile), h{r = 0) ~ lOOe, 
thus the maximal solution h{r,t — 2000 s), as a continuum construct, is at the limit of 
physical relevance. We also note that the experimentally observed spreading dynamics 
of a similar 5CB droplet typically stop after an hour, with R of the order of tens of 
millimeters (which corresponds to a pancake of volume V ~ 10~^ mm'^ and thickness 
30 nm). This behavior is due to short-range interactions which promote dewetting, 
and can not be captured by the TFE (fTTI) unless the expression |(9]) is refined. However, 
we do expect the TFE to capture the relevant properties of the contact fine region as 
observed experimentally by Cazabat et al |15[ [T6] , provided that the 30 nm thickness 
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Figure 3. Overview of the results of the numerical resolution of the thin film 
equation (TFE) with an "elastic" contribution to the disjoining pressure, i.e., a 
diffusive second-order term: a) three consecutive shapes adopted by a droplet of 
volume V = 2.36 x 10~^ mm'^ during a numerical spreading process (the oblique 
crosses indicate the points at which the solution h (r, t) falls below e = 30 nm, 
which is the characteristic thickness at the edge of the physical mesoscopic 
precursor [16)): b) same experiment: time dependence of the base radius R of the 
drop, as estimated directly from the tangent at the inflection point (see figure 1) 
and also, tentatively, from the volume conservation law V = ^R^O; c) comparative 
analysis of the estimates of R and 9 for the same cross-section snapshots as in 
(a): the tangent at the infiection point yields a sensible value for R, but 9 does 
not correspond to an apex-fitting spherical cap, even at t = Is. 
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plays the role of a cutoff, located in the asymptotic region of h{r, t) as is still the case 
for t = 2000 s on figures Hi and [3b. 

4-2. Review of the observables 

Figure [Sh makes it clear that the crossover from Tanner's law to a faster regime (at 
times of the order of a minute to an hour) coincides with a gradually lesser separation 
between the macroscopic and mesoscopic scales. In other terms, the extent of the 
mesoscopic "foot" of the drop as compared to the macroscopic "cap" is such that 
the apparent contact line is ill-defined. Another, somewhat unexpected fact shown by 
figure [Sfc is that measuring the contact angle 0inflcction(i) at the (mesoscopic) infiection 
point of h{r, t) does not give a good estimate of 9 for the (macroscopic) apex-fitting 
spherical cap, even in the Tanner phase (at t = 1 s) and despite a seemingly good 
separation of the scales. This observation, however, does not challenge our study, 
provided that we complement R, 0, -Rinflcction and ^inflection with additional observables 
that are consistent with those measured by Cazabat et al |15[ 116). 

As a matter of fact, the optical measurements in |15[ [T6] do not infer R and 
from the inflection point, which would have required a thorough reconstruction of the 
proflle at the contact fine for each snapshot. It was more expedient to track the edge 
of the spreading drop (located at an approximate thickness of 30 nm; we shall note 
this radius i?30nm) or the second interference fringe of the ordinary-extraordinary 
coincidence pattern (at a thickness of about 300 nm; we shall note this thickness 
-Rsoonm)- As for the contact angle, the slope at the contact line was averaged over 
the flrst 25 interfringes of the same pattern (between 300 nm and 5.3 pm). The 
numerical counterparts to these observable characteristics are presented in flgure 
Sbi, on a representative snapshot of the numerical spreading process. We note that 
the interference pattern used here is specific to nematics: the fringes correspond to 
coincidence between the ordinary and extraordinary rays, and the interfringe is about 
11 times larger than for the normal equal-thickness fringes. 

4-3. Time origin of spreading processes 

Both physical and numerical spreading processes can be seen as subject to initial 
conditions such as the deposit of a drop. Typically, shortly after a sufficiently 
compact deposit, fiow patterns appear at the edge of the drop and propagate 
throughout the initially static droplet, establishing Tanner's regime; at later times, the 
spreading crosses over to, e.g., a diffusive phase (for nematic droplets with antagonistic 
anchoring, as considered in this paper; c.f. figure [Sja). 

The physically relevant deposit is closely related to the subtle issue of choosing 
a time origin {t — 0). On one hand, the exact history of the deposit has no effect at 
the scale of the long spreading process. On the other hand, the power-law behavior 
typically observed in spreading is best represented in log-log diagrams (figure 2, 3b, 
3c, 4b, 4c), which are quite sensitive to the origin of t at small spreading times. 

Thanks to the robust presence of a Tanner stage at early spreading times, the 
dilemma is customarily resolved by describing the spreading in terms of the time 
elapsed since the effective origin of the Tanner phase: practically, for a sufficiently 
compact deposit, the origin of t is sHghtly adjusted so that, e.g., R{t) is well-fit by 
a power law at early spreading times. This may seem arbitrary but is in fact 
fundamental in the sense that for increasingly compact deposits of a given volume V 
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the spreading processes converge towards a well-defined limiting process, which (at 
least in the framework of thin film dynamics) precisely corresponds to a backwards 
extrapolation of Tanner's law. This convention is adopted for both the physical and 
numerical spreading processes presented in this paper. 



4--4- Scaling 

A prominent feature of the model TFE ifTTj) is that the equation can be scaled in terms 
of h, r and t (which corresponds to three degrees of freedom on cr, 77, D and V). If we 
know a function /io('^, t) that is a solution of 



dtho = - — -dr \ \hlrdr 
Vor [3 



-dr {rdrho] 
r 



Do-dr {rdrho) (15) 
r 



bearing the volume Vq = 271 h {r, t) rdr, then we can define 

k = {V/Vof"" (16) 



/ aor]D 
f OorjD^ 



1/8 



(17) 
(18) 



\crrioD^J 

and obtain a similar function h{r,t) — k^m^ho{k~^mr, k^^nt) which is a solution of 
l|lip with volume V. 

In the work presented here, besides the obvious fitting in terms of the volume V 

1 /8 

via k = {V/Vo) , we assumed that the surface tension a and the elastic coefficient 
K6^ were not significantly diflFerent from the values cr = 30 x 10^'^ N • m^^ and 
KS^ = 12 X 10~^^ J • m~^. We allowed, however, for an adjustment in terms of 
the effective viscosity 77, whereby m — 1 and n = r]o/r]. Indeed, the rheology of a 
nematic film with antagonistic anchoring conditions is not as trivial as the Poiseuille 
fiow in our model TFE: the effective viscosity must be intermediate between that of 
fiow-aligned 5CB molecules (30 x 10"'^ Pa • s) and that of fiow-orthogonal molecules 
(100 X lO"-^ Pa-s) [SniIlT]. This also affects the value of the effective diffusion coefficient 
D = The results of the fit are 77 = 70.5 x lO^^Pa-s, D = 1.136 x IQ-i^m^-s-i 
and V = 1.236 x IQ-^mm^. 

From the scaling factors k and n we conclude that the "characteristic time" of 

/ \ 1/4 

the crossover for the TFE ifTTj) scales as T = f ^V^/D^\ . For the numerical values 

yielded by the fit we have T ~ 8 x iCs (about 22 hours), which exceeds by far the 
duration of the physical experiment (2 hours). From the intersection of the asymptotes 
on figure 3b (at t ~ 700 s) and the values k = 1.23 and n = 0.425 we can extract a 
more quantitatively relevant time T — 5.7 x 10^ s, i.e., a couple of hours rather than a 
day. Both values are consistent with the fact that the crossover to a diffusive spreading 
process is far from complete at the end of the observation in [l6j . 



5. Comparison of physical and numerical spreading processes 

We shall now perform a quantitative matching between the physical and numerical 
spreading processes. We adjusted the scaHng of the numerical solution to 
accommodate the physical experiment, shifting the volume to V — 1.236 x 10^^ mm'^ 
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Figure 4. Quantitative fit of a physical spreading experiment (cyanobiphenyl 
5CB on silicon wafer) in terms of the numerical spreading process presented in 
figure 3: the volume V and viscosity r} act as scaling parameters of the maximal 
solution obtained previously (see text). Figure (a) illustrates the optically 
measured radii i?30nm (at the edge of the mesoscopic precursor) and ijgoonm 
(corresponding to the second of a series of interference fringes), as well as the 
optically measured estimate of the contact angle favcragod (averaged over 25 

fringes). Clearly ^inflection < ^?300nm < J?30nm and ^inflection > ^averaged- 

Figure (b) illustrates the good agreement for i?30nm and i?300nm- Figure (c) 
compares estimates of the contact angle, with a satisfactory agreement at late 
spreading times, and a discrepancy at early spreading times. 



and the viscosity to 77 = 7.05 x 10^^ Pa • s. Then, for this rescaled numerical 
experiment, we measured the same quantities as observed optically, namely 7?30nm, 
-Rsoonm and ^averaged (scc figure UK) . The results are presented in figures andHh. 

The best agreement was obtained for -Rsonm, over the whole range where the 
two experiments overlap. As for the measurements based on interference patterns, 
the agreement is less consistent. On one hand (figure Ilia), the radius i?3oonm of the 
second interference fringe agrees well with its numerical estimate. On the other hand 
(figure 3c) , at small spreading times the numerical ^average is in excess of the optically 
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observed slope averaged over the first 25 interfringes. The latter discrepancy may be 
due either to the low resolution of the fringes at early spreading times or to our failure 
to capture complex shear-thinning effects in the framework of the TFM. 

Although the agreement is quite satisfactory for the i?30nm observable, we must 
note that the numerical estimate of i?30nm is at the limit of applicability of the TFM. 
In our model we simply cut off the maximal solution of the TFE at the thickness /i = 30 
nm, and the dynamics of this cutoff line may differ from the actual dynamics of the 
mesoscopic precursor near the microscopic contact line. It would be more conclusive 
if the physical experiment had systematically provided the more robust observables 

-^inflection and ^inflection- 

6. Conclusion 

We have attempted an explanation of the abnormal spreading properties observed for 
small droplets of 5CB nematic liquid crystals |15[ll6j in the framework of the thin film 
model. This approach enabled us to illustrate both qualitatively and quantitatively 
the key trends in the spreading of nematic droplets: 

• the development of a large "foot" (mesoscopic precursor) , whereby the macroscopic 
and mesoscopic length scales are no longer well-separated; 

• the transition towards a faster spreading regime - determined by the antagonistic 
anchoring of the nematic at the interfaces - in which the thickness profile is 
essentially governed by a diffusion equation. 

The "acceleration", initially observed optically by Cazabat et al [HI [T^, was 
reproduced in a numerical spreading process (figure 3), which was used to fit the optical 
data (figure 2). The primary optical observable being the edge of the mesoscopic 
precursor, the agreement is satisfactory (figure 4b). 

We note that this post-Tanner regime is a priori not specific to nematic droplets. 
Similar crossovers to a faster spreading law than R ^ may be observed for regular 
liquids dominated by van der Waals forces, although perhaps not as readily as in the 
present case, where the spreading is driven by nematic elasticity. The fundamental 
result is that, for long-range substrate interactions, the droplet essentially becomes a 
diffusive film in the sense of Derjaguin [13] at late spreading times, and Tanner's law 
is gradually replaced with another law, determined by the substrate interaction rather 
than by capillarity. 

The work presented is a necessary complement to the quasistationary energetic 
approach as presented in [17], where acceleration is interpreted in terms of a dynamic, 
negative line tension r attributed to the apparent contact line (see figure 1). The 
latter framework is applicable if the macroscopic and mesoscopic length scales remain 
well-separated, i.e., if the bulk of the droplet is well-approximated by a spherical cap, 
and if both the vertical and lateral size of the mesoscopic region remain negligible. In 
this case, it is possible to isolate a line contribution, which resides in the mesoscopic 
"foot", yet contributes to spreading dynamics at the macroscopic scale. However, as the 
droplet spreads, it eventually adopts a characteristic bell shape, and the capillarity- 
dominated spherical cap ceases to be a good approximation. At even later stages, the 
droplet may reach the state of a mesoscopic pancake, which can not be resolved by 
the macroscopic model at all, unless the dimension is lowered to a planar geometry. 

By contrast, the thin film model appears to provide a more robust description 
of complete wetting situations, especially in the late stages of spreading. In order 
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to account for the emergence of pancakes, we aim to provide the thin film equation 
with suitable boundary conditions that would account for the phenomenology of the 
microscopic contact line. As a prospect of future work, we also note that our current 
implementation of the thin film model does not accurately describe the non-Newtonian 
rheology of antagonistically anchored nematics. In future studies we may refine the 
notion of effective viscosity and allow for a more accurate modeling of the spreading 
dynamics. 
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